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A coupling between the spacetime geometry and a scalar field involving the Euler four-form can 
have important consequences in General Relativity. The coupling is a four-dimensional version of the 
Jackiw-Teitelboim action, in which a scalar couples to the Euler two-form in two dimensions. In this 
case the first order formalism, in which the vierbein (or the metric) and the spin connection (or the 
affine connection) are varied independently, is not equivalent to the second order one, in which the 
geometry is completely determined by the metric. This is because the torsion postulate = is not 
valid and one cannot algebraically solve the spin connection from its own field equation. The direct 
consequence of this obstruction is that the torsion becomes a new source for the metric curvature, and 
even if the scalar field is very slowly varying over cosmic scales as to have no observable astronomical 
effects at the galactic scale, it has important dynamical effects that can give rise to a cosmological 
evolution radically different from the standard FRWL model. 

1 Introduction 

A century after Einstein's formulation of General Relativity (GR) no one doubts that gravitation is a 
manifestation of the geometry of spacetime: matter and energy curves spacetime on which matter and 
energy move. It is by looking at the trajectories of particles and light that one can find out the spacetime 
geometry and thereby infer the local energy-matter distribution. Similarly, by studying the shape and 
evolution of the universe one infers the distribution of energy and matter at the cosmic scale. Thus, at 
the galactic scale, the velocity distribution of stars suggests the presence of (dark) matter of an unknown 
nature; at the cosmic scale, the accelerated expansion of the universe points to the existence of a form of 
vacuum (dark) energy with exotic properties {p < — l/3p). 

These puzzling results have landed over the past decades on top of long standing riddles like the 
unexpected smallness of the cosmological constant and the quantum understanding of gravity. Even if 
these conclusions were not the result of inaccurate measurements, or wrong assessments of the data, it 
might be a healthy attitude to consider variations on the assumptions of the theoretical framework and 
to explore the scenarios resulting from such variations. 

1.1 The standard cosmological model 

The current view is that the universe is roughly described by the standard Friedmann-Robertson-Walker- 
Lemaitre model, also referred to as big-bang cosmology. The main ingredients in this model are: 

• Symmetry (principle of equivalence): Spacetime is a differentiable manifold endowed with a bundle of 
tangent spaces and the laws of physics are invariant under local Lorcntz transformations [U El [3] . 

• Dynamics: Extremizing the Einstein-Hilbert action 



with respect to the metric g^j^, yields Einstein's equations (EE), 




(1) 



(2) 
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• Cosmological principle: The spacetime metric is a solution of EE that admits a global slicing into 
homogeneous and isotropic spacelike three-surfaces. 

The validity of this model is based on astronomical observations which seem consistent with the 
following facts: 13.7 x 10^ years ago, the universe was hot, small, fairly homogeneous and isotropic; the 
spatial section of the universe is currently quite flat and homogeneous at large scales; the universe is 
very near the critical density; the expansion of the universe seems to be accelerating. The theoretical 
challenge is to explain these observations with a consistent dynamical model. 

In order to model dark energy and/or dark matter the simplest option may be to include some exotic 
fluids that contribute to the right hand side (RHS) of the EE (see, e.g., [4]). An alternative could be to 
change the dynamics of geometry, defined by left hand side (LHS) of EE, replacing the Einstein-Hilbert 
Lagrangian by some scalar density function of the metric, the curvature and the torsion f(g,R,T). For 
arbitrary /, this in general brings in new degrees of freedom for the graviton through higher order 
derivatives (see for instance [5l[6j[7]), but an equally serious concern is the fact that there seems to be 
no principle from which the function / can be derived. This problem is like that of the cosmological 
constant but not just for one arbitrary parameter but for infinitely many. 

In the explorations described above a clear distinction between what goes into LHS (geometry) and 
RHS (matter-energy content) of ^ is implicitly assumed. This scheme is a direct extension of the experi- 
ence in mechanics and electrodynamics, where the sources (RHS) determine the evolution of the relevant 
dynamical variables -positions and field configurations- (LHS). This point of view is so ingrained in our 
way of thinking about spacetime that we infer the matter content of the universe from the motion of par- 
ticles. If the geodesies are straight lines, we assume spacetime to be flat and there is no matter anywhere. 
Conversely, if the geodesies are not straight lines, we infer that there must be matter somewhere, even 
if invisible. This, however, is not necessarily true. It has been shown that there exist forms of matter 
coupled to gravity that do not necessarily curve spacetime [8]. 

Here we show that the converse may also fail to be true: there could be sources for local curvature 
that are not necessarily forms of matter, the source of curvature might be geometry itself, and this may 
result in a radically different cosmological evolution. As we shall see next, this situation results from a 
minimal modification of GR, showing how sensitive on the assumptions of GR the resulting cosmology 
the can be. 

1.2 Second order formalism and torsion 

Another important simplifying assumption postulated by Einstein in the standard approach to GR is 
the vanishing of torsion. According to this condition, the affine connection is symmetric and defined 
by the Christoffel symbols, so the metric is the only fundamental field to be varied in the action. The 
Einstein-Hilbert action is a functional of the metric and its variation yields second order equations for 

On the other hand, in the first order formalism, the action is varied with respect to the vielbein 
e" and the Lorentz connection oj""^ independently, in a similar scheme as that proposed by Palatini [9]. 
In this case torsion is not assumed to be zero but is determined from the field equations. In four- 
dimensional matter-free gravity the two formalisms are equivalent because in the first order approach, 
the field equations still imply T° = 0. So, on shell the torsion-free condition is satisfied in both cases. 

Moreover, in the 4D and matter-free case, the equivalence is also true off-shell because the equation 
T'^ = 0, obtained by varying the action with respect to w"'', can be algebraically solved for a; as a function 
of e° . By the implicit function(al) theorem, this means that the reduced action obtained by substituting 
Lu — a;(e), the second order action I[e] = I[e,ui{e)] is completely equivalent to the original first order one, 
even away from the classical extrema [10] . Thus, in four-dimensional matter- free GR, the torsion- free 
condition does not restrict the dynamics around classical solutions (see, e.g., pTl fT2] l. 

If T" does not vanish, its presence would be revealed by fermions, as they would follow different 
geodesies from the those for spin-0 or spin-1 particles, an effect that may not be significant for the 
current experimental accuracy (see, e.g., [13 )• However, even though postulating = might seem 
harmless and could pass the experimental tests in four dimensions, it is not a logical necessity and it is 
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theoretically unsatisfactory to impose this constraint that is not respected in the presence of ordinary 
matter, like spin 1/2 fields. 

Another, perhaps more important aspect of torsion, is that it may act as a source for the metric 
curvature. This may have important cosmological consequences. An example of this effect is presented 
here. 



1.3 The scalar-Euler form coupling 

The modification we consider here consists of adding a coupling term between a scalar field and the Euler 
density, 

^eabcdR'"" A W\ (3) 

where i?"'' = dw^^'+w^cAw^'' is the curvature two-form for the Lorentz connection uj°'b- For constant <j) the 
modification ([s]) is irrelevant, as it adds a constant to the action proportional to the Euler characteristic of 
spacetime. So one can expect that if is slowly varying over cosmological distances, the Einstein equations 
should not be affected at planetary scales. However, even if at present the effects of the modification could 
be negligible, the additional term changes the field equations by adding a new equation, EabcdR""^ ^R'^'^ = 0, 
and modifying the other equations in a form that could significantly affect the evolution near the big-bang 
as well as at the final stages of the universe's expansion. 

The scalar field (j) could be an effective field whose origin is not relevant for the present discussion. 
With the possible exception of the Brout-Englert-Higgs particle, all scalars in nature are composites of 
spin-1/2 matter fields, whose origin could be found in the underlying fundamental microscopic theory. 
Mathematically, the scalar coupling to the Euler form can be found in the dimensional compactification 
of a higher-dimensional gravitational theory of the Lovelock family, where the first nontrivial correction 
to the Einstein-Hilbert action in dimensions D > 5 is the Gauss-Bonnet term. For instance, in five 
dimensions this term is 

eABCDpR^'' A i?*^^ A e^. (4) 



Reducing to four dimensions, and using 



(5) 



produces ^ |14| along with other non- minimal couplings involving torsion explicitly jl5| . 

The scalar-Euler density coupling ([3]) is the natural extension to four dimensions of the two-dimensional 
Jackiw-Teitelboim gravitation theory [16l I17j. whose gravitational Lagrangian is the product of a scalar 
field and the two-dimensional Euler form, cfieabR"'''- 

This coupling could also be seen as a particular choice of the most general extension of scalar-tensor 
theories originally proposed by Horndeski [T5] , recently revived in the context of Galileons \W[ [501 121] > 
where the combination - R^.^ + R''^^"' RaPt^u) is one of the possible extensions of gravity 

that yields second order field equations for the metric (the comparison between using V{(j)) and simply (p 
is discussed in the last section). 

The main difference between our work, originally discussed in [22], and the Horndeski approach is 
that we adopt the first order formalism where e° and w"'' are independent fields independently varied 
in the action. This makes a great difference because in the presence of scalar-Euler form coupling the 
first and second order formalisms are no longer equivalent. The contribution of ([3| to the variation of 
the action with respect to o;"^ implies that torsion does not vanish if is not constant. Moreover, this 
equation is not algebraically solvable for the connection as a function of e° and (j). Hence, there is no 
local expression uj°'^{e, (f>) to substitute in the action. Thus, the presence of the term ^ prevents writing 
the action only in terms of the vierbein (or the metric) and (j). 

On the other hand, in the second order form, the field equation eabcd.R"'^ ^R'^'^ — 0, obtained by varying 
(f), implies the vanishing of the Kretschmann invariant (i?"^^'^i?a/9/^i/)- This equation represents a strong 
constraint on the geometry, incompatible with most classical solutions of GR, such as the Kerr black 
hole, which rules out a coupling like ([s]) in the second order formalism, unless additional modifications 
to the theory are introduced. This is not the case in the first order formalism since the condition ^ 
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establishes a relation between purely metric part of the curvature and torsion which, as we shall see in 
section 3, it does not necessarily eliminate interesting options. 



2 First order action and field equations 

Let us consider the first order form of the Einstein-Hilbert action with cosniological constant in four- 
dimensions, including the scalar-Euler density coupling ([3||^ We will also allow the possible presence 
of other forms of non geometric matter represented by some unspecified fluid, whose dynamics is not 
relevant for the present discussion. The action reads 

/[e, w, 0] = ^ J ea6cd(i?°'e=e'' - ^e'^e'e'e'' + ^R^'R'^) + J (6) 

where A is the cosmological constant, and L„i represents the non geometric matter. Varying with respect 
to e", a;"'' and 0, yields the following field equations 

Se" : eabcd{R'"'-^e-e'')e' = 16nGT^ (7) 

Soj''' : eabcd{T''e'' + dcj)R-'')^0 (8) 
Scj) : eabcdR''''R'''' = 0. (9) 

where is the stress-energy 3- fornj^ defined by SL^ — Se"^ A t^. Eq. ([t]) has the same form as the 
Einstein equations with the important difference that since now ([8| implies that torsion can no longer be 
assumed to vanish, i?"** is not the usual Riemann curvature of standard GR. This is best seen splitting 
the Lorentz connection as w"'' = (D"'' -I- where w"'' is the torsion-free part, 

de"^ + cS-, A e" = ^ = eS(V^K), (10) 



and the contorsion k""^ determines the torsion, = A e^. In (10 1, is the inverse vielbein and V 
is the covariant derivative in the Christoffel connection. In this scheme, the Lorentz curvature splits into 
de usual torsion-free part of Riemannian geometry, R"''', and the torsion-dependent terms, 

^ab ^ ^ab ^ Q^ab _^ ^ ^c6^ (11) 

where D is the covariant derivative in the torsionless conection a)"''. In the presence of torsion the K-terms 
become sources for the Riemann curvature i?"'' in the Einstein equations. Using (11), Eq.(|7| becomes 



^abcd ^ 



R'' " ^e'e^) - -eabcd(D«;^^ + AtV^-^^)e'' + 16ttGt, , (12) 



where the torsion terms on the RHS add to the other matter sources represented by r^. 

In the metric formulation, the field equations for the same action would have been completely different. 
The metric expression corresponding to (|3| is 



i?2 „ 4i?^-i?^, + R'^'^'fR^^^p ) y=^, (13) 



which has been noted to give rise to second order field equations [181 [21]. However, varying with respect 
to g^jy under the assumption of vanishing torsion yields, instead of ([T]) and ([8|, the Einstein equations 
plus quadratic curvature corrections and terms containing both first and second order derivatives of (j). 
It can be easily seen that the resulting sets of equations of both formulations are inequivalent. 



^Henceforth wedge products of forms will be assumed. 



^The stress-energy 3-form is related to the stress-energy tensor as *T£j = Tcds'^, where r'^d = |e| ^e^ '^^^^j.d"'^ 
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3 Cosmological consequences 



We now consider a spacetime foliated by a family of isotropic and homogeneous three-dimensional spatial 
slices (Bianchi type I, V & IX), as described by the standard Big-Bang cosmology, and we examine the 
solutions of the field equations ([7]-|9| in this case. The metric. 



1 — kr^ 



'{d0^ +sm^ 9dip^) 



(14) 



admits six global Killing vectors {£^giiu = 0), associated to three spatial translations, ^(i), and three 
rotations, wher^ 

= - kr'^d^^ , and = Xid^ . - Xjd^.. (15) 

We assume torsion and the scalar field to have the same isometries as the background spacetime, 
£iT°' iiu = and £^0 — 0. It is straightforward to see that these conditions imply (j) = 4>{t), and 
the only nonvanishing components of T"' f^, are 



T^'ev = 2/(t)a(i)rVl - kr^ smO 
2f{t)a{t) sine 



rriU 



Vl - fcr2 



Tr 
rt 



'2.f{t)a{t) 
\/l — kr'^ sin ( 



et 



rpip 



hit) 



(16) 



where h(t) and f{t) are functions of time to be determined by the field equations. 

As shown in Appendix A, substituting the expressions for the metric an torsion into the first order 
field equations ([7]-[9|, yields 



k 2 A 

3 3 
k 



2{U + HU) + + ^ - r - K = -n^p, 



-f){U + HU)~2fU{f + Hf) = 



(17) 
(18) 
(19) 
(20) 
(21) 



where H ~ d/a is the Hubble function, — SttG, pm and are the energy density and pressure 
described by the matter Lagrangian, and U = H + h. Taking the covariant derivative of ([T]) the stress 
3-form is found to satisfy the on-shell condition Dr^; = T'^m^^Tc or equivalently. 



K ide Tc 



(22) 



were D = d + cj denotes the covariant derivative in the torsion-free connection. As seen in t he appendix, 
here = h{t)e^ and k^'^ = — /(i)e^"^if e^, which implies that the right hand side of (22) vanishes 



identically. This implies the continuity equation for the matter content of the universe, p„ 
Pm) = 0, which is ultimately related to the fact that the metric ( [l4| and the torsion (161 satisfy Bianchi 
identities TlR"-^ 



dR^b + uj^cR^b - R^cOJ^b = 0, and DT° 



Equations (17) and (18) correspond to the modified Friedmann equations and can be cast in the 
standard form: + ^ ~ — ifP ^^'^ 2ij + + ^ — A = ~K^p, where p = pm + Pt and p = p„i +pt 
combine the density and pressure of ordinary matter and the effect of torsion. Thus, the effective density 



''The usual relation between cartesian (xi, X2, x-j) and polar (r, 6,1/3) coor dinates is understood. 

■'This is not an obvious property of our ansatz, because the expression is only a consequence of the tensorial character 
of fj,jy in a space of that admits the Kihing vectors (llSlt and does not depend on being the torsion tensor. 
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and pressure produced by torsion can be identified as 

3 



Pt 
Pt 



1 



{2Hh + h'' - r), 

{2h + mh + h^ - f). 



(23) 
(24) 



This identification is made more compelling by the fact that px and pT can be shown to satisfy the 
corresponding continuity equation, pT + iH{pT + Pt) = 0. 

3.1 Solutions for cold pressureless matter 

We consider ordinary matter described by dust Pm — 0. The continuity equation p„i + 3iJpm = 0, implies 
Pm — Po(ao/a)'^)j where po represents the present density of the universe, and og is the present value of 
the scale factor which can be normalized to one. In order to solve the system, we first note that since 



equation (21) represents the Euler invariant [23 , it is locally a total derivative. Indeed, (21 1 can be 
written as 



dt 



3" 







whose integral is 



3 Va2 J ) 



(25) 
(26) 



where C is an integration constant. Assuming / 7^ in (20) implies </> = (2C/) ^, and replacing (19) in 
(18) allows integrating U{t) as 

Uit) = tanh [^/A/2it - to)] , (27) 
for A > 0, where to is an integration constant that sets the origin of time. For A < the solution is 



Uit) 



tan 



V-A/2(t - to) 



and for A = 



U{t) = {t- to)- 



(28) 



(29) 



Note that these results do not depend on po or k. Hence, from (17), (20) and (26) a, and / are obtained 



and h can be finally read from (19). The solutions are different for each sign of A: 
i)A > 



A 



a{t) 

h{t) 

fit) = VA 



cosh 



log 



1 
A 

1 /A 



3 V 2 



sinh 
coth 



A/2(t - to) 
^/K/2it - to) 
v/A72(t - to) 



3V2C 



An Po 



1/3 



sinh 



/A72(i - to) 



cosh 



fKj2it - to) 



3y2C- 


\/Kk^Po tanh"^ 


yA72(t - to) 


i^/2C- 


Ak^po tanh 


y/A/2it - to) 



fcsinh^/^ 


v/A72(t-to) 




cosh^ 


7172 (t - to) 


(3\/2C 


- V Ak^po tanh 


7X72 (t - to) 


y/3 



2C- 3Ctanh2 


yA72(t - to) 


+ ^K^potanh^ 


yA72(t - to) 


6C - y2AK2potanh 


v/A72(t-to) 





1/2 



(30) 
(31) 
(32) 



(33) 



ii)A < 

ait) 

hit) 



: COS 



V-A/2(i - to) 



3%/2C 

y=A72(t - to) 



V^k'^Po 

y=A72(t - to) 



1/3 



1 



log 



A 

1 /-A 



sm 



v/-A/2(< - to) 



cot 



V-A/2(t-to) 





\/—Ak'^Po tan"^ 


v/-A/2(t-to) 


3\/2C + 


\/— Ak^Po tan 


V-A/2(t - <o) 



fit) = V^l- 



k sec^ 


v/-A/2(t - to) 


sin2/3 


y-A/2(t - 


-to) 


(3^C + \/^K2potan 


V-A/2(t - to) 


y/3 



2C + 3Ctan 


y-A/2(t - to) 




g^K^po tan^ 


v/-A/2(t-to) 


6C + V-2Ak2pq tan 


V-A/2(t-to) 





1/2 



(34) 
(35) 
(36) 



(37) 



iii)A = 

ait) 

hit) 
fit) 



= ^(t-to)'/'(«:Vo-3C7(t-to))i/3 
= lit - to)^ + <l>0 

K^PO 



3(t-to)(K2po-3C(t-to)) 

( 22/3fc 



K^po - 9C(t - to) 



1/2 



(t-to)^/3(K%-3C(t-to))2/3 3(t-to)2(K2po-3C(t-to)) 



(38) 
(39) 
(40) 

(41) 



The integration constant to corresponds to the origin of the time coordinate, which may be taken as the 
end of the radiation-dominated era. For A > 0, tg is chosen as zero. For A < 0, a convenient choice is 
one for which a(0) = 0, that is to — — 7r/\/— 2A for A < 0, and to = for A = 0. 

In Fig. ^ we observe that for C 7^ 0, for A > the scale factor decreases from an infinite value, 

collapsing at a finite time if P > 1; in the limiting case P = 1, the scale factor tends asymptotically to a 

finite constant. For A < there is an initial singularity, if A < there is an infinite expansion at a finite 

,2 

cosmic time, while with A = there is always a Big-Crunch at tc — ^^3^; note that, counterintuitively, 
the time of collapse tc grows with the density. 

The special case C = must be treated separately (the sign of a is irrelevant since only appears 
in the metric), 
a) For A > 0: 

( cosh^ 



• a 

• h ^ tanh 

• P ^ fccosh 
b) For A < 0: 



v/^A72(t-to)]) 



1/3 



V-A/2(t - to) 

4/3 



v/-A/2(t - to) 
1/3 



tanh 



v/-A/2(t - to) 



a - {^cos^ yA72(t-to) ) 



(a) for A>0 (b) for A<0 

a(0 a(t) 




t 



(c) for A=0 

a(t) 

1 r - Q=o.l 

Q=0.5 




Figure 1: Scale factor for A positive, negative and zero respectively. For A ^ 0, i and a{t) arc measured 
in units of |A|~^/^, for A = a{t) is measured in the same units of t. The graphics show the behavior of 

a(t) for different values of the addimcnsional parameter P = '^/^^^^ for the case A ^ 0, or for A = 

2 

the parameter Q = ^^^fP- measured in the same units of t 



(a) for A>0 (b) For A<0 




t 



Figure 2: Function /^(t) for A positive, negative an zero respectively. For A 7^ 0, t is measured in imits of 
|A|~^/^, while P{t) is measured in units of |A|, for A = f^{t) is measured in units oft~^. The graphics 
show the behavior of /^(t) for different values of the spatial curvature k and for particular values of the 
parameters po. C, and A. Continuous line corresponds to the case fc = 1, dashed lines is for k = 0, and for 
k = —1 we use dotdashed line. Note that, in general, there are some values of t for which f^{t) becomes 
negative. 
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• h ^ tan 

c) For A = 0: 



A/2(t-to) 

A/2{t - to) 

2/3 



tan^ 



/A72(t-<o) 



4 Discussion and summary 

A minor variation of the form of coupling discussed here would be to change (jj — > V {((>). This only means 
that the function that couples to the Euler density is not a fundamental field, and offers some flexibility 
for the field equations because in the equations it replaces d(j) by V'{4>)d(p. This results in different 
scenarios for V'{(t)o) = and V'{4>o) 0, where (f)o is a classical configuration. We have not discussed this 
case in detail here as it does not add a very significant conceptual difference to the previous discussion, 
although it may entail important cosmological consequences that deserve careful analysis. The specific 
nature of the field 4> is not very important, its origin could be exotic or mundane, but the fact that is 
present makes a huge difference. 

One could try to refine the model by assuming a dynamical scalar field </>. For instance, assuming 
matter to be described by the Klein-Gordon Lagrangian L„i = ^(dcj))'^, the density and pressure of the 
scalar field become p = p = -^(p- The modified equations are completely integrable (see Apendix B for 
details). However, a numerical analysis beyond the scope of this work would be needed in order to obtain 
detailed information about the solutions. 

Although our model might be regarded as not very realistic, it shows how a slight modification in the 
form of the action, having negligible effects in a small region (e.g., at the scale of the solar system at the 
present age of the universe), can give rise to dramatic changes in the cosmic evolution. Let us summarize 
the main features of this very simple scenario: 

• For A > 0, the solution does not have an initial, point-like singularity. The universe starts from an 
(infinitely large scale, contracting to a minimum size before starting an era of accelerated expansion. 

• For A < 0, there is an initial point-like singularity as in the standard FRWL cosmology. At late 
times, the universe also experiences an accelerated expansion, reaching an infinitely large scale factor in 
a finite cosmic time, t = 7r/\/— 2A. 

• For A = the universe also starts with an initial singularity, but in this case there is no accelerated 
expansion, the universe, instead, contracts and collapses in a finite time, tc = ^^3^^, that is proportional 
to the present density of the universe . 

• The time evolution near the big bang is highly sensitive to the sign of A. The solution for A = 
cannot be obtained by taking the limit of A — > from the A 7^ solutions, however its behavior is very 
similar to the case A < near the big bang. 



• In this model, torsion was introduced through a scalar field non-minimally coupled to gravity, plus 
the adoption of the first order formalism, a slight modification of the standard assumptions of General 
Relativity. Thanks to torsion, the universe undergoes an accelerated expansion without requiring either 
dark energy or negative pressure fluids. In fact, the nonvanishing torsion contributes to both "matter 
density" (23) and "pressure" (24). In this model, the universe can undergo an accelerated expansion for 
both signs of the cosmological constant, while in standard GR with A < 0, an accelerated expansion can 
only take place in the presence of exotic matter, or other exceptional conditions |24| . 
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Apart from the Euler form, there exist other topological invariants like the Pontryagin and the Nieh- 
Yan[321 Unj forms that could also be analogously coupled to a scalar field. We are aware of several 
explorations in this direction, in the second order formalism |261 1271 128] . Membranes of co-dimension two 
or more can also couple to Chern-Simons forms [29. ,30] , bringing in torsion as well. 

As seen in Fig. ([2]) for A > and small t — to, as well as for A < and large t — to, p becomes 
negative. This situation may seem catastrophic, but since / corresponds to the completely antisymmetric 
part of the contorsion, it docs not affect the geodesic equation of classical particles and does not couple 
to spin zero or Gauge vectors fields. This component of the torsion only affects spin ^ particles through 
the coupling in the Dirac Lagrangian |13) 

/:, = -^T„^^V^r"'5^^ = ^/Vir°r5v.. (42) 

This interaction term couples torsion to the chiral current and is added to the torsion-free Dirac La- 
grangian Ld — |(V'r''D^?/' — Dp'!/)r^'0) — irmpip. If / becomes imaginary, the corresponding loss of 
unitarity and violation of the current conservation can be interpreted as giving rise to particle creation. 

Among all the possible solutions discussed above, the case A < and k — seems closest to the present 
cosmological evidence: a big-bang singularity with an accelerated expansion and a negligible presence 
of torsion effects at late times. The current limited understanding of the evolution of the universe at 
large scales urges a revision of the standard assumptions of cosmology. The case presented here is just 
one example of how dramatically different are the resulting scenarios if certain logically unnecessary 
assumptions we take for granted -like the absence of torsion or the metric formulation of GR-, are 
removed. Another example where this kind of assumptions are removed is presented in [31] 



Appendix A 

With the convenient choice of the vierbein e" 
the torsion 2-form become^ 



dt and = yf^^^dr, = a{t)rd9, = a{t)r sm6d(p, 



rpQ 





h{t)e'e" + f{t)e'jKe'e''. 



(43) 
(44) 



From the equation = K^be'' we identify the contorsion one-form to be = h{t)e^ and k^"' — 
—f{t)e^'^KC^- On the other hand, the nonvanishing components of the torsionless part of the spin con- 
nection are Cj'^^ — He^ , cj^^ = — (ar)^-'^e^\/l — kr'^^ u)^^ = —{ar)~^e^y/l — kr"^ and Co^^ = — (ar)^^e'^ cot 9] 
where H = a/a is the Hubble function, and therefore the spin connection reads 



U! 



01 



\/l — kr'^ 



13 



UJ 



23 



Vl — kr^ 

( 

ar 

cote 



(45) 
(46) 

(47) 
(48) 



The Lorentz curvature i?"** = do;"'' + ijj"-c^'^^ takes the form 



01 



= [{H + h)+ H{H + h)]e°e' + f{H + h)e^jKe-'e^ 



R 



ij 



^Here I, J, K = {1,2,3} are Lorentz indices. 



(49) 
(50) 
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The stress 3-Form is related to the stress tensor by = ■^^labcT^de'^s^s'', and considering a perfect fluid 
r'd = diag{—p,p,p,p) we obtain 



To 

Tl 



(51) 
(52) 



Finally, replacing (|43) lii) |49) [50]) into the field equations 0-([9]), one obtains Eqs.(|17|)-(|2f 1). 



Appendix B 



One way to refine the model is to assume a dynamical scalar field <j). In that case, the matter Lagrangian 



^d^ * d^. The stress tensor is Tab — da 
- 4>{t), the density and 



hilabdc4>d''(j) (where we defined da — E^d, 



is Lm(0) 

and since (f> = (/>(t), the density and pressure of the scalar field are given by p = p 
B is also modified as eabcdR""^ R'^'^ = An^d * d(j). Assuming again that / 7^ 0, (26) is replaced by 



all 

^(p^- In this case, 



U 



(- 



-r - 



12U 



C_ 



(53) 



Substituting (17) in (18 191 again reduces the problem to quadratures that can be integrated for U{t). 
If |A| ^ K, then 



t-tn 



A2 - k2 



arctan 



2U 



. V-A - VA^ 

The cases A — ±k must be handled separately, giving 



arctan 



2U 



t-tn 



2U 



tanh ^ 


' 2U' 







4C/2 - K 

respectively. For A = the solution is 



& t-tn 



2U 



arctan 



v/-A + \/A2^^ 



2U 



+ K 



2K{t — to) — arctan 



1 - 




— arctan 


1 + 




1 , 

+ 2 log 















{Vsu- 



/kY + K 



iVsu 



Finally, combining this with Eqs. (I7|, (19), (53) and = the system can be completely integrated 
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